A magnetic tomography of a cavity state 
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A method to determine the state of a single quantized 
cavity mode is proposed. By adiabatic passage the quantum 
state of the field is transfered completely onto an internal Zee- 
man sub-manifold of an atom. Utilizing a method of Newton 
and Young (Ann. Phys. 49, 393 (1968)), we can determine 
this angular momentum state uniquely, by a finite number of 
magnetic dipole measurements with Stern-Gerlach analyzers. 
An example illustrates the influence of dissipation. 

PACS Nos. 42.50.-p, 42.50Vk 

The state of a quantum mechanical system is com- 
pletely specified by its density operator p. It is a fun- 
damental as well as an important practical question 
of quantum mechanics to devise measurement schemes 
which allow a complete determination of p. By a se- 
quence of repeated measurements on an ensemble of iden- 
tically prepared systems the state has to be character- 
ized operationally. In quantum optics this topic of com- 
plete state determination has received recently consider- 
able attention in the context of characterizing nonclassi- 
cal states of the radiation field, and states of atomic and 
molecular motion Jl| |ic| . 

In a seminal paper, Vogel and Risken [|lj have pointed 
out that the state of a single mode of the radiation field 
(equivalent to a one-dimensional harmonic oscillator) can 
be found by a tomographic techniques and correspond- 
ing experiments have been performed by Raymer and 
coworkers |2| . The central idea of quantum state tomog- 
raphy is based on a reconstruction of the density ma- 
trix p from measured quadrature probabilities p x (9) = 
(x, 9\p\x, 0). Here \x,0) is a rotated eigenstate of the 
quadrature operator xg, i.e \x, 9) = 1Z{9)^ \x) = e iea a \x) 
with a and lowering and raising operators of the os- 
cillator. Alternative schemes have been discussed in the 
literature under the name of state endoscopy or by 
introducing discrete Wigner- functions R|. Very recently, 
ideas were developed for a state determination of ions 
moving in harmonic trapping potentials l^jic}]. 

In the present paper we discuss a new scheme to mea- 
sure the density matrix of the radiation field of a single 
quantized cavity mode by a magnetic tomography. It is 
based on combining ideas we have developed in the con- 
text of quantum state engineering of arbitrary Fock-state 
superpositions in a cavity by adiabatic passage |ll|l , with 
a tomography of atomic angular momentum states by 
Stern-Gerlach measurements originally proposed in fl^| . 

It is well know that the adiabatic change of a Hamil- 



tonian interaction transforms initial energy-eigenstates 
into eigenstates of the final Hamilton operator |13[ . This 
method to map quantum states is applied in various con- 
texts of molecular- and atomic physics |Il],|l4],[l5| . It is 
also the key mechanism for transferring the state of a 
quantized cavity mode onto the internal state manifold 
of an atom. 
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FIG. 1. Evolution of a degenerate two level atom (angu- 
lar momentum: J 9 — > J e = J g — 1) coupled to a quantized 
cro-polarized cavity mode. First, it passes through the profile a 
classical beam (cr+) and then, with a delay, through the cavity. 

According to Fig. |l|, an atom passes adiabatically 
|]IT| , ^5| through the spatial profile of a classical de- 
polarized laser beam [Rabi-frequency: f2(f)] and, with a 
spatio-temporal displacement r > , through the profile 
of a quantized, 7r-polarized cavity mode [atom-cavity cou- 
pling: g(t — r)]. We assume that the electronic structure 
of the atoms corresponds to an optical J g — > J e — J g — 1 
dipole transition. 

The coupled atom-cavity system evolves according to 
the time-dependent Hamiltonian 
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tion frequency of the cavity mode, respectively. In terms 
of atomic basis states and Clebsch-Gordan coefficients 
C a lm'm e , the atomic de-excitation operators A ,Ai are 



defined by 
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The state space spanned by this Hamiltonian has the re- 
markable feature that it can be decomposed into invari- 
ant sub-spaces H = 0^1_ 2 j W l '. Due to angular mo- 
mentum conservation, it is only possible to couple angu- 
lar momentum states to a finite number of photon states 
by means of a unitary evolution (Eq. [I]) . 

U v = {\m g \ <J g ,0<n = v + Jg + m g \\ 
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One element of the sub-space TL V is of particular interest, 
i.e. a linear combination involving only ground states 
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By an appropriate choice of coefficients a v m ,i.e. 
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this normalized state (a 1 ', 



TV) becomes also an 
eigenvector of the Hamilton operator in the interaction 
representation (derived from Eq. (|l])) and has a zero 



eigenvalue. According to the adiabatic theorem 13 1, this 
eigenvector approaches a stationary eigenstate of the cor- 
responding Schrodinger-equation, if the time dependent 
change of the Hamilton operator during the total interac- 
tion time T is much less than the characteristic transition 
frequencies (O max T, \/N max g max T > 1). Furthermore, 
if the delay and shape of the pulse sequences are chosen 
such that 
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then this is a mapping process that only permutes states 
up to a sign change s v — sigo.(a v _ J _ v (+oo)). 
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In other words, a coupled atom-cavity density operator 
p\AC) ^ na ^ can factorized initially into a pure atomic 
state and a field state containing less than 2 J g photons 
will be mapped to a product of atomic ground state su- 
perpositions and the cavity vacuum 
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With reverse adiabatic passage, an internal atomic state 
is prepared uniquely by reading out the cavity state. 

The complete characterization of such an angular mo- 
mentum state by a number of magnetic dipole measure- 
ments was described in Ref. [jl2). It requires to detect a 
set of physical observables that are proportional to 

{|m| < J,\s\ < 2J\\P m (6,ip s )}. 

Here P m (6,ip s ) represents the projector onto a rotated 
state \J,m,TZ(6,(p s )e z ) = D^ JS> (TZ)\J,m,e z ), s enumer- 
ates an arbitrary set of 4J+1 azimuthal angles tp s and 
is a constant inclination. This method can be imple- 
mented, for example, by the unitary evolution of an an- 
gular momentum state in a homogeneous magnetic field 
B(8, <p s ) = \B\ n(9, ip s ) oriented differently, each time the 
measurement is performed, and by using a conventional 
Stern-Gerlach analyzer Fig. |^. 



O 

• 1— I 

o 

e 



O 

o 

u 
Oh 




P„ t (9.<P) 



FIG. 2. Setup of a Stern-Gerlach experiment with an addi- 
tional homogeneous magnetic field B{0,ip) = \B\n(6,ip) induc- 
ing spin precession around the axis n(6,ip) 
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From the J-dimensional representation of the rotation 
group D^ J \lZ) or the Wigner matrices S J \Q) ||i7| ; one 
finds 



\J,m,H{0,(p)e z ) = e- mip d^(6)\J, ni e z 



n=—J 



Hence, the occupation probabilities are given by 



(6) 



For simplicity, it is assumed that the cavity and the ex- 
ternal laser have a common frequency u c — ujl and are 
detuned from the atomic resonance by A = w eg -Wi. The 
resulting atomic density operator can be determined ei- 
ther by solving the master-equation (Eq. (|l0|)) or alterna- 
tively, by averaging over a number of simulated quantum 
trajectories |is| ]. 
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This linear equation relates density matrix elements 
((2 J + 1)(2J + 1) real numbers) to measured prob- 
abilities that are positive numbers. By determining 
(2 J + 1)(4J + 1) different probabilities, this seemingly 
over determined set of linear equations has a unique so- 
lution that is positive definite. For \l\ < J , < w < J — I, 
one finds 
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In case of an equally spaced array of azimuthal angles 
ip s , i.e, — 7r < ip s = s ^jt^ < \s\ < 2J, the quantity 
X{9) is the discrete Fourier transform of the measured 
probability tableau 
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In contrast to systems with continuous degrees of free- 
dom, the reconstruction algorithm of Eq. (H) is faithful 

if inclination angles 6 are avoided where d|/o(#) vanishes 
(i.e. the zeros of an associated Legendre polynomial). 
Most detrimental to this state tomography is the loss of 
cavity photons during the adiabatic interaction. In con- 
trast, spontaneous atomic decay is of minor importance 
as the adiabatic eigenstate is formed by a ground state 
superposition. To examine the influence of dissipation, 
we have coupled the atom-cavity system to an environ- 
ment 11 1 and obtained the following master equation for 




the density operator p. 



dt P 
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where T and k denote the spontaneous decay rate and 
the inverse cavity life time, respectively. 

In the interaction picture representation (derived from 
Eq. |l|), the effective, non-hermitian Hamiltonian, intro- 
duced above, is given by 



FIG. 3. a) Real part of the initial cavity density matrix p^l 
vs. photon number n,n; b) atomic ground state occupation 
A a pA\ + KCi c pal, (10) probabilities p m (9,ip s ) vs. magnetic quantum number m and 
azimuthal phase index s of an atom that moved adiabatically 
through the cavity; c) atomic density matrix p^f— vs. magnetic 
quantum number m,m, inverted tomographically from the data 
represented in Fig. Bb (see text for parameters). 



In Fig. the results of the mapping and reconstruction 
process are shown for an initially pure cavity state 
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The real part of the initial cavity density matrix, i.e. 
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p( c ^(to) = \^(ta))c{'ip(to)\c vs. photon number n,n is 
shown in Fig. |3|a. In order to map this cavity state 
onto a Zeeman submanifold, we assumed a sufficiently 
large degeneracy (J g = 4 — > J e = 3). Both fields 
are tuned to the atomic resonance A = 0. All fre- 
quencies are scaled to the spontaneous decay rate of 
the atomic excited state T, as the peak Rabi- frequency 
Sl max = 50 r and the maximal cavity coupling constant 
.9max = 30 r. The time dependent Gaussian turn-on 
(beam) profiles were of identical shape FWHM = 1 [r _1 ] 
and had a relative delay of t = 0.65 [r -1 ]. To com- 
plete the adiabatic passage, a total interaction interval 
(to = -1.07 [T- 1 ], t end = 1.72 [r- 1 ]) was chosen. The 
cavity decay rate was set to re = 0.01 T which implies a 
decay probability of Pc — K(N)(t en d ~ to) ~ 0.15 (for 
contemporary cavity QED experiments in the optical or 
microwave regime see Ref. |19[| ). 

From the final density matrix, the atomic ground state 
occupation probability (Eq. Q) was evaluated with ip s = 

S 47TT' N ^ 2J and 

Pm(^,<Ps) = (J,m,n(^,ip s )\ p^j (tc„d) \J,m,n(^,<p s )). 

The result is depicted in Fig. ^b. 

Subsequently, we applied the tomographic inversion 
(defined in Eq. |§|) to these data to obtain p^f^, (Fig. ||c). 
Direct comparison with the simulated density matrix 
Psim^cnd) shows that the inversion procedure induces no 
error. The additional features that appear in the mapped 
quantum state are of physical origin. The population 

(A) 

p\ 4 is solely due to the spontaneous emission of a cr_ 
photon, as this state is otherwise not coupled to the dy- 
namics at all. On the other hand, the satellite peaks 
that appear in vicinity of the original coherent superpo- 
sition state are caused by the decay of the cavity state 
during the adiabatic mapping process. However, the re- 
semblance with the original cavity state is striking. 

In summary, we have studied a method to map the 
state of a single quantized cavity mode adiabatically 
onto a finite dimensional degenerate Zeeman submani- 
fold of an atom that passes through the resonator. Sub- 
sequently, we characterize this state by a number of re- 
peated Stern-Gerlach measurements on identically pre- 
pared atoms as outlined in Ref. By a full quantum 
mechanical calculation, including spontaneous emission 
and cavity decay, we have shown that this method yields 
a faithfull image of the original, a priori unknown cavity 
state. This method is not limited to the measurement of 
pure states but may be applied also in case of statistical 
mixtures. 
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